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Abstract: A new lattice estimate of K ^ 2tt transitions claims, contrary to all 
other computations, that the hadronic matrix element for the gluon penguin operator 
Qq has opposite sign and, in addition, is much larger than the vacuum saturation 
approximation. We comment under what conditions (if any) it is possible to reconcile 
this lattice result with the experimental value of e'/ e . The dramatic impact of new 
physics in the kaon system that seems to be required is not easily accommodated 
within our present theoretical understanding. 



Keywords: Kaon Physics, CP violation. Lattice QCD. 



1. A new lattice estimate of K ^ 27i transitions using domain- wall fermions 
claims [1] , contrary to all other computations, that the hadronic matrix element for 
the gluon penguin operator Qq has opposite sign and, in addition, is much larger than 
its vacuum saturation approximation. This surprising result comes about because 
of the contribution of the so-called eye-contraction diagrams. The resulting value 
for the CP violating parameter e'/ e is of the opposite sign and almost one order of 
magnitude bigger than the current experimental determination: 

= (2.1 ±0.46) X 10-^ (1) 

which is obtained by averaging over the KTeV [2] and NA48 [3] preliminary results 
as well as the older 1992-93 experiments (E731 [4] and NA31 [5]). 

If this lattice result (taken at its face value) will stand further scrutiny it raises 
serious questions on our understanding of electroweak physics in the kaon system. 

In this letter we would like to discuss under what conditions (if any) it is possible 
to reconcile a large and positive {2 7r\QQ\K) with the experimental value oie' /e . We 
consider two possible scenarios: 

• A modification of short-distance physics that changes both the sign and the 
size of the Wilson coefficient of the gluon or electroweak penguin operators. 
This can in principle be achieved by tampering with the initial matchings of 
the various coefficients while preserving the standard basis of operators. As we 
shall discuss, the behavior of the rcnormalization group equations (RGE) for 
the relevant efi'ective lagrangian force us to rather extreme changes in order to 
achieve the desired effect. Even though supersymmetry provides a framework 
for potentially large effects, to construct a model in which such large devia- 
tions from standard physics are present in e' while conspiring to be invisible 
everywhere else involves a contrived set of assumptions. 

• An enlargement of the standard operator basis. We limit our discussion to the 
case of the chromomagnetic penguin operator because it does not affect the 
renormalization of the Wilson coefficients of the other dimension six operators 
(allowing us to draw model-independent conclusions) and it has been already 
shown to be a promising candidate for new-physics effects in e' /e . Given 
the current estimate of this operator's matrix element, short-distance changes 
alone — although potentially very large — are not sufficient in bringing e' /e in 
the experimental ball-park. A final assessment requires a lattice evaluation of 
the hadronic matrix element of the chromomagnetic operator. 

Both scenarios call for a dramatic impact of new physics in the kaon system and they 
are not easily accommodated within our present theoretical understanding. More 
exotic extensions of the standard operator basis and more extreme scenarios require 
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a detailed model-dependent discussion which is beyond the scope of the present 
letter. ^ 



2. Let us fix our notation by introducing the A^* = 1 efi^ective four-quark 
lagrangian for mc < < 

= - ^) CMQ.iix) - + r C,{ij)Qi{ii)} , (2) 

i=l,2 i=l,10 

where Ag = VqdV*g and r = — At/A^. In the discussion that follows we will assume 
as the standard model (SM) reference value for the CP phase ImA^ = 10"^. The 
renormalization group Wilson coefficient Cj(//) are known to the next-to-leading order 
in Us and ae [7, 8]. 

The standard basis of effective operators is 

Q2 = (sm)v-A Mv-A ' Q2 = (sc)v-A (cC^)v-A ' 

Q3,5 = (srf)v_A (^?)vtA ' 
Q4,6 = {Sadl3)y_j^T,q{q(5(l 

Q7,9 = 1 (so?)v-A Eg (gg)v±A > 

Qs.lO = I (^arf/3)y_^X^qeq(^^go,)v±A, 

where V ± A stands for 7^(1 ± 75) and iq is the value of the electric charge of the 
quark q = u,d, s, c. 



3. Taking the hadronic matrix elements of the operators in eq. (3) obtained by 
the lattice calculation [1], and given their standard model Wilson coefficients, e'/e is 
completely dominated by the contribution of Qq and is therefore large and negative 
(in dramatic disagreement with the experiment). 

In order to reconcile the lattice and the experimental result, we must somehow 
compensate for this dominant contribution. The ratio e'/e is determined by the sum 
of isospin / = and 1 = 2 amplitudes. Let us examine in turns possible modifications 
on these two classes of contributions. 

A first possibility, in the / = amplitude, is that the Wilson coefficient of the 
Qq, at the scale at which the matrix element is computed (about 2 GeV), is changed 
with respect to its standard model value in order to reproduce the experimental 
result. Assuming that new physics only modifies the initial conditions (at niw) of 
the RGE, we need to know how these changes are propagated by the RGE down to 

^Possible effects of dissipation and loss of coherence in the kaon system relevant to s'/s have 
been discussed in ref. [6]. 
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Figure 1: The Wilson coefficient Cq{2 GeV) as a function of the parameters r2 and rg 
defined in eq. (4). 

the hadronic scale. This has been discussed for the whole operator basis in section 
VII of ref.[9]. 

In order to keep the discussion as model independent as possible, we parameterize 
all deviations from the SM matching conditions Ci{mw) in terms of the parameters 

r, = a{mw)/Cr{mw). (4) 

In general, a Wilson coefficient at the low scale receive contributions from both a 
multiplicative and additive renormalization, the latter arising from QCD-induced 
operator mixing. In particular, the Cq{^) is dominated in the standard basis by the 
additive renormalization induced by the mixing of with the Q2 operator. 

We have plotted in Fig. 1 the Wilson coefficient Cq{2 GeV) as a function of the 
parameters r2 and r^. The dark (red) band represents the values of Cg for which the 
experimental result for e' / e would be recovered (leaving all other Wilson coefficients 
unchanged), that is 

Cq{2 GeV) 
Cl^'{2 GeV) 

As it can seen by inspection, this possibility is realized either by a drastic reduction 
of C2{Tnw) (the left side of Fig. 1) or a large enhancement of CQ{mw) (the right 
side of Fig. 1). However, the value of 6*2(2 GeV) is severely constrained by the 
CP-conserving amplitude Aq{K 27r), which in the same lattice calculation is 
reproduced up to a factor of two. Moreover, it is difficult to substantially change 



-0.2 < < -O-l ■ (B) 
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IJ, — mw Ci C2 C3 C4 C5 Cq C-j Cs Cg Cio 



Tree ^/ 
Tree + g V V 
Tree + 7 



V V V V 




Table 1: Contributions to the one- loop matching of the AS = 1 Wilson coefficients at 
= niyy. Non- vanishing contributions to Cg and Cio arise via the QCD renormalization of 
the operators and Qg, respectively. The label "Tree" stands for tree- level W exchange, 
^g,7,z for gluon, photon and Z-induced penguins, while B for the t^-induced box diagrams. 

C2{mw) by means of new physics because this Wilson coefficient is due to tree level 
1^-exchange. 

Looking then at Fig. 1, and given that C2{mw) cannot differ too much from 
its standard-model value, it is possible to reproduce the experimental e'/e only by 
enhancing 0^(171]^) by more than a factor of twenty. However, such a large enhance- 
ment can hardly take place without affecting other processes and we shall come back 
to it after discussing the electroweak penguin. 

Leaving the / = contribution alone, we can still compensate for the large and 
negative result by acting on the the 1 — 2 contribution. Here the dominant operator 
is Qs- Given the lattice estimate of the matrix element of this operator, agreement 
with the experiments would require 



The RGE analysis shows that Cs{2 GeV) varies proportionally to C'j{mw) (recall 
that Cii{mw)=Q) and, as one can see from Table 1, C'j{mw) receives contributions 
from photon and Z penguins. 

To date, the best limits on the CP conserving component of the Z-penguin 
operator, ReVz, are provided by the decay [10, 11, 12], whose branching 

ratio is measured to be [13] 



even though the constraint on KeVz is not as accurate as the experimental preci- 
sion because of the theoretical long-distance uncertainties related to the two-photon 
component [14]. 

The cleanest constraint on \Vz\ comes from the decay tt^vv, which is 

currently measured to be [15] 




(6) 



B{Kl ii^ii-) = (7.18 ± 0.17) X 10 



-9 



(7) 



B{K+ 



TT+Pu) = 4.2 tli X 10 



-10 



(8) 
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Due to the lack of further evidence, the branching ratio in eq. (8) is going to be 
reduced by a factor of two or three [16]. 

Taking the standard model expectation B{K+ n+uu) = (0.8 ±0.3) x lO'^" [17] 
as a reference point, a numerical analysis shows that we can at most modify ImPf^ 
by a factor of sixteen. 

Barring the unlikely possibility of independent and widely different effects in the 
effective {sd)v-A vertex of the photon and Z penguins, eq. (8) rules out the enhance- 
ment required by eq. (6). Analogous considerations hold for the other electroweak 
operators which would entail even larger deviations of their Wilson coefficients from 
the standard model values. By the same token, also the enhancement by a factor 
of twenty of the gluonic penguin coefficient CQ{mw) — which we argued is necessary 
in the 7 = amplitude in order to reproduce the experimental value of e'/e — is 
difficult to accommodate. 

A notable exception to the above argument is the presence of gluino-induced 
flavour-changing neutral currents in supersymmetric models. The gluino induced 
AS = 1 transitions are suppressed in Z-penguin diagrams compared to gluon or 
photon penguins by a factor of 0{m\/rn\) (for a detailed discussion of this effect 
which follows from gauge invariance see ref. [18]). As a consequence, the bounds on 
Z-penguins are not effective on gluino-induced gluon and photon penguins. However, 
the gluino contributions to the standard = 1 penguin operators are indirectly 
constrained by the sharp bounds on the gluino-induced /S.S — 2 box diagrams relevant 
to - K° transitions [19]. 

Gluino-box diagrams can also play a direct role in AS — 1 amplitudes. In 
ref. [20] , it is shown that gluino-box diagrams may induce a potentially large isospin- 
breaking contribution to the electroweak penguin Qg while satisfying all other bounds. 
However, the large factor required by eq. (6) implies a rather large mass splitting 
between the right-handed squark isospin doublet (m^^ — mj^ ~ 1 TeV) together 
with a large gluino mass in order to evade the bounds from Atuk and e. 



4. An operator not usually included in the standard model analysis of e'/e is 
the chromomagnetic penguin 

Qii^^s [TTidR + rUsL] a^, G^" d , (9) 

where R{L) = (1 ± 75)/2 and G^^ is the gluon field. 

The matrix element {2'k\Qii\K) has been computed in the chiral quark model 
and shown to arise only at 0{p^) in the chiral expansion and to be further suppressed 
by a rri^/m^ factor with respect to the naive expectation [21]. For these reasons, 
even though the Qn Wilson coefficient receives a large additive renormalization, its 
standard model contribution to e'/e has been shown to be neghgible [21]. 
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This is no longer true if the Cabibbo-Kobayashi-Maskawa (CKM) suppression of 
its CP violating component can be lifted without violating other bounds. A clever 
example of it is discussed in ref. [22], where the standard model factor ImAt ~ 
0(sin^ do) is replaced in a supersymmetric framework by a CP violating squark mix- 
ing of 0{sin.9c)-i thus introducing a potential enhancement of the chromomagnetic 
Wilson coefficient by three orders of magnitude. This solution allows for a contri- 
bution to e' /e at the 10^^ level by saturating the known bounds coming from CP 
violating phenomenology (discussions of the implications of these bounds on various 
supersymmetric models can be found in ref. [23]). 

This enhancement of Cu is still not enough to compensate for the huge neg- 
ative Qq contribution to e' /e of the lattice result; in fact, keeping fixed the gluon 
penguin contribution, agreement with experiment would require a matrix element 
for Qii larger by about a factor of ten. Actually, the leading O(p^) chiral quark 
model estimate of the hadronic matrix element may receive potentially large 0{p^) 
contributions if the accidental ml/mj^ suppression is replaced by mf^/A^. To further 
assess this possibility, it would therefore be interesting to have an estimate of the Qu 
matrix element from the same lattice approach that produces the large and positive 
matrix element for the Qq. 

Clearly, extensions of the standard model can also involve new effective operators 
beyond the standard basis of eq. (3) and Qu. However, in this case very little can be 
said without a complete re-analysis of e'/ e . Similarly, scenarios in which the CKM 
matrix is taken to be real and CP violation arises only in the new-physics sector can 
in principle be invoked but again require a detailed model-dependent analysis before 
being considered a viable alternative. 

5. In conclusion, while a combination of the above scenarios may make some 
of the requirements less severe, a viable model which avoids the phenomenological 
bounds discussed is forced to rely on a contrived choice of parameters. It is neverthe- 
less remarkable how supersymmetry can provide a framework for potentially large 
effects on e'/e while satisfying all present data. 

It is also fair to add that the lattice result by means of domain-wall fermions 
in ref. [1] must stand further scrutiny and corroboration before concluding that the 
standard model is facing its most dramatic challenge to date. 

We thank F. Borzumati and M. Neubert for useful discussions and comments. 
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